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ABSTRACT 

Does star formation proceed in the same way in large spirals such as the Milky Way and in smaller chemically younger galaxies? 
Earlier work suggests a more rapid transformation of H2 into stars in these objects but ( 1 ) a doubt remains about the validity of the H2 
mass estimates and (2) there is currently no explanation for why star formation should be more efficient. M 33, a local group spiral 
with a mass ~ 10% and a metallicity half that of the Galaxy, represents a first step towards the metal poor Dwarf Galaxies. We have 
searched for molecular clouds in the outer disk of M 33 and present here a set of detections of both '^CO and '^^CO, including the only 
detections (for both lines) beyond the R25 radius in a subsolar metallicity galaxy. The spatial resolution enables mass estimates for the 
clouds and thus a measure of the A'(H2)//co ratio, which in turn enables a more reliable calculation of the H2 mass. Our estimate for 
the outer disk of M 33 is A'(H2)//co(i-o) ~ 5 x 10^° cm"-/(K km s"') with an estimated uncertainty of a factor < 2. While the '-"'CO 
line ratios do not provide a reliable measure of A'(H2)//co> the values we find are slightly greater than Galactic and corroborate a 
somewhat higher A'(H2)//co value. Comparing the CO observations with other tracers of the interstellar medium, no reliable means 
of predicting where CO would be detected was identified. In particular, CO detections were often not directly on local HI or FIR or 
Htt peaks, although generally in regions with FIR emission and high HI column density. The results presented here provide support 
for the quicker transformation of H2 into stars in M 33 than in large local universe spirals. 
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1. Introduction 



I Several recent pape rs dLerov et al] 120061 : iGardan et al.l 120071 : 
. iGratier et alj|2010bh have suggested that the rate of transforma- 
' tion of molecular gas (H2) into stars is higher in small chemi- 
[ cally young galaxies - those with low metallicities. The latter 
- articles devoted much attention to whether the conversion factor 
' from CO to H2, A^(H2)/Ico, usually expressed in H2 molecules 
, cm"^/(K km s"'), could be severely underestimated. In this arti- 
' cle, we present sensitive measurements of the CO emission from 
] the outer disk of M 33, allowing us to estimate Virial masses 
,; for isolated molecular clouds far from the center as well as some 
' '^CO line measurements, providing a further check on the physi- 
cal conditions in the Giant Molecular Clouds (GMCs) in M 33. A 
further goal is naturally to understand the mechanisms of molec- 
ular cloud formation and the outer disk provides a means to ex- 
plore physical conditions unexplored by earlier work which only 
presented molecular cloud data in the inner disk where the stellar 
mass surface density dominates that of the gas. 

In particular, we present data on a series of mid to outer disk 
clouds in M 33, including an interarm GMC with no measured 
associated star formation (Ha or Far-IR emission) and a GMC 
beyond the R25 radius. Because the stellar population in M 33 is 
rather young, the R25 radius (30.8' i Paturel et al. 200J ), defined 
in B band, corresponds to an extremely low stellar surface mass 
density, well below that of the gas at the same radius. In M 33, 



Galaxies: Local Group - Galaxies: evolution - Galaxies: ISM - ISM: Clouds - Stars: 



the radius at which stellar and gaseous surface densities are equal 
is about 0.5 R25 whereas in large spirals the stars dominate until 
about the R25 radius. 

It has long been known that the amount of CO emis- 
sion per unit star formation or per unit H2 is lower, some- 
times much lower , in galaxies with subsolar metallicity (e.g. 
iRubio et aP Il99lh . In the solar neighborhood, the conver- 
sion factor A^(H2)/Ico is of order 2 x 10^"cm"^/(Kkm s"') 
with a strong increa se from the center to the outer disk (e.g. 
'Sodroski et al.|[l995h . This radial variation is not only true for 
the Galaxy but also for other spirals where it has been studied 
dBraine et al.ll997allbh . The amount of CO emission per H2 mass 
varies not only with metallicity but also with the radiation field, 
which can act to raise the CO emission by warming the gas. If 
the metallicity is low the effect is opposite because the lower 
dust content and lesser degree of self-shielding by molecules 
other than H2 reduce the size of the CO-emitting regions with 
respect to the H2 cloud size, such that the CO emission princi- 
pally traces the central regions of molecular clouds. Because the 
CO lines are optically thick, moderate changes in metallicity and 
radiation field with respect to Galactic conditions (solar vicinity 
or molecular ring) can be expected to only generate moderate 
changes in the A^(H2)/Ico ratio. 
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Fig. 1. lOjum image of M 33 ( iTabatabaei et al.ll2007h using a logarithmic scale to show the outermost features. The R25 radius is 
illustrated with a yellow ellipse. The inner black ellipse indicates half the R25 radius, where gaseous and stellar surface densities 
are roughly equal. The positions observed in '^CO are shown in yellow in the outer parts and in black in the inner parts. Positions 
observed in '''CO are shown in red except for the two innermost points (M33_l 1 and M33_12) which are shown in blue in order to 
be visible. All the positions ob served in '^CO ha ve also been measured in '^CO. The white contour shows the 1.5 x 10^' cm"^ HI 
contour at 17" resolution from lGratier et alJ (l2010a) . 
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The standard techniques to estimate the H2 mass, and thus 
A^(H2)/Ico ratio, are through the dust emission and an assumed 
dust-to-gas mass ratio, the "virial" mass from the cloud size and 
Hne width, and from optically thin tracers such as '^CO. We dis- 
cuss the latter two methods in Sections 4. 1 and 4.2. The resulting 
N(ii2)/Ico is then used to estimate the star formation efficiency 
(SFE), defined as the ratio of the star formation rate per unit H2, 
inM33. 

M 33 has a metallic it y a factor 2 be low solar 
dRosolowskv & SimonI 120081: iMagrini et al] |2009|) . a mod- 
erate UV radiation field, and is a spiral galaxy despite having 
a mass ~ 10 times lower than the Galaxy. It thus represents a 
first step towards the still smaller, lower metallicity, irregular 
galaxies in the Local Group and beyond. Based on the above, we 
assume an Oxygen abundance of \2+log(0 1 H) - ^A-Q.03Riij,c. 
The distance to M 33 is assumed to be 840 kpc and we adopt 
an inclination angle of 56° and a position angle of 22.5° as in 
iGardan et all (l2007h . 



2. Observations 

The observa t ions presented here are a follow-up to the 
iGardan et"an (l2007h mapping of a large part of M33. All data 
were taken with the 30meter telescope run by the Institut de 
RadioAstronomie Millimetrique (IRAM) on Pico Veleta near 
Granada, Spain. Observing runs were in April, August, and Nov. 
of 2006. 

The i-CO and '^^CO /=1^0 and 7 = 2^1 transi- 
tions were observed with the "AB" receivers. These receivers 
are dual polarization and the beam-splitter enables simultaneous 
observation of the two transitions. Data were taken under good 
conditions, with system temperatures (T*) of 300 - 400 K ex- 
cept for the '''CO(l-O) where the average system temperature 
was ~ 200K. As a result, data reduction was simple: bad chan- 
nels were eliminated, a zero-order baseline (constant value) was 
subtracted, and spectra were averaged for each position. 

The positions observed are shown as yellow and black 
circles in Fig. 1 for the '^CO and red or blue circles for the 
'^CO. The size illustrates the IRAM CO(l-O) half-power 
beamsize. Generally, both the 7 = 1 — > and 7 = 2 — » 1 
transitions were observed simultaneously. The data are pre- 
sented using the main beam temperature scale. The main beam 
efficiencies are estimated to be 0.72, 0.74, 0.53, and 0.54 for 
respectively the lines at 115, 110, 230, and 220 GHz (see 
http://www.iram.es/IRAMES/mainWiki/Iram30mEfliiciencies) 
and forward efficiencies 0.95, 0.95, 0.91, and 0.91. 

Table 1 provides the positions at which we observed '"CO. 
Tables 2 and 3 give the '^CO 7=1^0 and 7 = 2 — > 1 fluxes 
at each detected position and the nns noise level at IMHz res- 
olution for the undetected positions. The results of fitting gaus- 
sian line profiles to the spectra, in order to robustly measure line 
widths and velocities, can be found in Table 4. Table 5 gives the 
'^CO fluxes and '-^'^CO line ratios for both transitions. 



3. Molecular clouds in the outer disk of M 33 

In this article we present CO and even '^CO detections beyond 
the R25 radius of M 33. These are the first such detections for 
a sub-solar metallicity galaxy. Only a few detections have been 
made in the outer parts of galaxies: in the Mi lky Way, using the 
rotation curve to estimate the dist ance (e.g. iBrand et alJ 1987t 
iDigel et al.| [ T994l). in NGC 4414 ( iBraine & Herpinll20o4 . and 
in NGC 6946 dBraine et alJl2007h . One of the issues is whether 



large quantities of molecular gas could have escaped detection 
in the outer disk (iPfenniger et al.|[l994l) . In very nearby galaxies 
such as M33, molecular clouds can be spatially resolved, allow- 
ing the use of the Virial theorem and isotopic lines to (roughly) 
assess their masses and resemblance to inner disk clouds. This 
is particularly interesting for small, chemically young, galaxies 
Hke M 33. 

After the map ping of the northern part of M 33 by 
IGardan et al.l (l2007h . we wanted to integrate more deeply at po- 
sitions where we thought CO might be found, with the goal of 
testing whether we could predict from other means if molecular 
gas were present. As molecular clouds form out of the atomic 
gas, which dominates in the outer disk, one of the criteria was 
the presence of high HI column density. Dust emission is fre- 
quently linked to CO emission but in the outer parts of M 33 
the FIR and MIR emission was often too weak to be detected by 
Spitzer (or earlier satellites), despite the substantial HI column 
densities. A series of positions was observed to a lower noise 
level, and in CO (1 -0), rather than in the CO(2-l) mapping by 
IGardan et al.l (l2007h : the coordinates are given in Tab. 1 and they 
are shown in Fig. 1 as circles. Figures 2, 3, and 4 present the 
CO(l-O) and HI spectra at these positions on the 70fj.m image 
with HI column density contours. The HI data is at 17" reso- 
lution fro m a mosaic of VLA C and D configuration fields de- 
scribed in lGratier et al.l (l2010ah . These outer positions were cho- 
sen because either star formation was present from the Ha im- 
age (lGreenawalllll998l: iHoopes & WalterbosllIOOOl) or high HI 
antenna temperatures (either narrow lines or high t otal co lumn 
densities) were found in the iDeul & van der Hulstl ( Il987h data. 
The FIR data was not available at the time of the observations 
yet the spectra show that the CO detections are not always at the 
HI or FIR/Ha maxima although they tend to be part of larger 
HI structures with Ha emission. At higher sensitivity, the FIR 
emission would probably be detected as well. Here we show the 
70jum data but the 160jum images do not show more extended 
emission. Figures 2, 3, and 4 show the sources labelled M 33 18, 
20, 22, 23, 24, 26, and 27 - all of which are norfli of Dec. 3 1 :05. 
Only deep NIR images detect the stars at these distances from 
the center 

Fig. 5 shows the CO(l-O) and CO(2-l) spectra of the outer 
disk clouds M33.2, M33_3, M33^, and the CO(l-O) spectra of 
the two detected positions in M33_20. Fig. 6 shows the 8 clouds 
for which '^CO emission was detected as well. Among the '^CO 
detections is cloud M33_18, beyond the R25 radius (see Fig. 3 
as well). This is the first such detection and '''CO is detected 
in both transitions. It is also of interest that the line widths de- 
crease significantly from the inner regions (bottom of figure) to 
the outermost parts (top). The narrow lines show that we are only 
observing a single cloud, allowing us to try to use the Virial the- 
ore m in the next section to estimate the cloud masses. Cloud 4 
(the IGardan et al.l 120071 "Lonely Cloud") is very well-detected 
but is in an interarm region with little or no detected star for- 
mation despite the strong CO emission. Clouds 4 and 18 have 
been observed with the Plateau de Bure Interferometer and will 
be discussed at length in a dedicated article. 

While there is a general link between CO emission and HI 
and FIR {e.g. lOfim) emission, the CO emission of clouds 4 
and 18 appears surprising even with the benefit of hindsight. 
The 70//m emission, for example, is only very marginally de- 
tected at the positions of these strong CO detections. However, 
something must provoke the condensation of atomic gas clouds 
to form H2. In the case of cloud 18, the HI spectra are 
complic ated so perhaps the H2 is the result of merging HI 
clouds (iBrouillet et al.l 119921: iBallesteros-Paredes et all 119991: 
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Fig. 4. CO(l-O) spectra (yellow) and HI spectra (red dotted) 
of source M33_18 and major offsets on 70/^m image with HI 
contours at 0.8, 1.2 (dashed), and 1.6 xlO^' cm"^ (dotted). The 
wedge shows the lOfiin image scale in MJy sr"'. The CO maxi- 
mum is actually at the (-10.5,0) position shown in Fig. 6. 



iHennebelle & Peraulll 119991: iHeitsch et al] l2005h . In the other 
molecular clouds in the outer regions of M 33, this seems less 
likely because the HI spectra are not wider or more complicated 
than in many places where CO is not found. 
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Fig. 5. Spectra of the 4 clouds where only '^CO emission was 
detected. The 7 = 1 — » transition is in black and the 7 = 2 — » 1 
is in red (dashed). The spectra of M33_20 show the offset (150,0) 
with a solid line and the dashed line is the (120,0) offset, both in 
CO(l-O). The ,Y-axis shows the velocity in km s"' and the y-axis 
the antenna temperature on the main beam scale. The order is 
in increasing declination from bottom to top, corresponding to 
increasing galactocentric distance. 



4. Determining the I-I2 mass 

4.1. via the Virial theorem applied to individual clouds 

Historically, cloud masses have been estimated assuming that 
molecular clouds are bound under their own gravity. Galactic 
GMC s show a relationship between size and linewidth dLarsonl 
|1981|) which is difficult to explain in other ways. The masses 
calculated are usually referred to as "Virial masses", twice 
as large as marginally gravitationally bound clouds. In addi- 
tion to GMCs, d iffuse clouds are also observed in the Galaxy 
(iPolket al.lll988h and these are probably unbound. The scales 
we observe in M 33 correspond to GMCs, justifying the use 
of the Virial method. Since the only force opposing gravity in 
the standard calculations is the internal motion, measured by 
the line width, magnetic support could increase masses further if 
clouds really were stable long-lived "Virialized" objects. Since 
the idea of virialized clouds came out, it has become clear that 
GMCs are not long-live d entities b ut rather have "lifetimes" 
of order 10^ years (e.g. iKawamuraet al. 2009). However, the 
"Virial" masses have generally been found to be consistent with 
other means of measuring G MC masses, such as dust measure- 
me nts (|Sodroskietal.l 119951) or gamma-ray ob servations (see 
e.g. Combeslll99ll) . Using the expression from iDickman et aP 

mm . 



Mv 



AV-D 

2yG 



where AV is the half-power velocity width and D the cloud 
diameter ( y/AA/n for other shapes). The y parameter is prin- 
cipally related to the broadening of the '^C O line due to its 
high optical depth (see iDickman et al.|[l986l for details). Using 
this, the cloud masses (Myir) are given in Table 4, along with 
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Fig. 6. Spectra of the 8 clouds where '-^CO emission was de- 
tected. The 7 = 1 — > transition is in black and the 7 = 2 — > 1 
is in red. The '^CO antenna temperatures are multiplied by 10. 
The dotted red line is the '''CO(2-l) and the dashed black line is 
the '^CO(l-O) spectrum. The ,Y-axis shows the velocity in km^"' 
and the y-axis the antenna temperature on the main beam scale. 
The order is in increasing declination from bottom to top, corre- 
sponding to increasing galactocentric distance. C loud 8 is most 
likely a blend of two clouds at different velocities (lGardanl2007h 
and in the spectra above there is a clear difference between the 
CO(l-O) line shape and center compared with the (2-1) transi- 
tion. 



and masses calculated using a 
2/(K km s-i) factor (Miun,)- 



cloud velocities, line widths 

^(H2)//cO(l-0) =2x10 cm "/l,rs. KJii s ' ) lacLOi (.Wilum 

The Virial masses given are upper limits in the sense that we 
have assumed that the clouds are as big as the beam, 85 pc in 
diameter, whereas this is likely an overestimate of the size of a 
single cloud (e.g. Solomon et a l. 1987 ). Furthermore, multiple 
clouds would broaden the CO lines and again cause an overes- 
timate of the mass. However, the 7 = 2 value may be overes- 



timated and is used for consistency with iBraine et al.l (Il997ah . 
These authors also found, through basic radiative transfer calcu- 
lations and for modestly subsolar metallicites, that a low radia- 
tion field (as for the clouds in the outer disk of M 33) could coun- 
teract the deeper CO photo-dissociation expected with decreas- 
ing metallicity, such that the CO-emitting part of a cloud was not 
necessarily smaller in an outer disk environment. Nonetheless, 
the virial masses given here are higher (possibly due to the large 
assumed size) in most cases than galactic GMCs of similar line 
width (using Mvii ~ 365AV'^Mq, eq. 6a from ISolomon et al.l 
Il987h . Attributing an uncertainty is difficult; the sources of er- 
ror are chiefly the size, the y parameter, and the hypothesis of 
"Virial" equilibrium (gravity alone). The size and y likely lead 
to an uncertainty of a factor 2 at most. The presence of CO- 
bright isolated clouds in the outer disk is evidence that they are 
robust, and thus gravitationally bound, structures, comforting the 
Virial hypothesis. If the clouds we identify were simply due to 
the fluctuating turbulent field then we would expect to find many 
weak lines for a single stron g one (if indeed stron g lines could 
be formed this way) and the iGardan et al.l ( l2007h observations 
suggest this is not the case. The relatively strong '■'CO emission 
shows that the '^CO line is highly optically thick, also arguing 
for a non-transient nature. The spectra of some clouds, particu- 
larly in the inner disk, are complicated and single-gaussian fits 
lead to significant uncertainties in the linewidth and thus Virial 
mass. These are a minority, however, so we estimate an uncer- 
tainty of at most a factor 2 on the value given below. 

With the exception of the two bright inner disk clouds, 
M33_ll and M33_12, which are consistent with a galac- 
tic A^(H2)//co(i-0) factor, and the very uncertain sources, the 
virial masses are consistent with a A^(H2)//co(i-0) ~ 5 x 
10^°cm"'^/(Kkm s"') factor (the average for clouds 4, 5, 7, 9, 
10, and 18). 
Table 1. Positions observed in CO. 



Name 


RA 


Dec 


SRA 

(") 


6Dec 

(") 


M33.1 


01:34:17.67 


31:05:12.0 


0.0 


0.0 


M33_2 


01:34:06.13 


31:03:09.7 


0.0 

-20.5 

-0.0 

20.5 

-0.0 


0.0 

-0.0 

-20.5 

-0.0 

20.5 


M33.3 


01:34:34.86 


31:00:42.9 


0.0 

-20.5 

0.0 


0.0 

0.0 

-20.5 


M33.4 


01:34:16.81 


30:59:30.0 


-10.5 

-10.5 

-10.5 

0.0 

0.0 

0.0 

10.5 

10.5 

10.5 

-21.0 

21.0 

0.0 

-0.0 


-10.5 

-0.0 

10.5 

-10.5 

0.0 

10.5 

-10.5 

0.0 

10.5 

-0.0 

0.0 

-21.0 

21.0 


M33.5 


01:34:03.06 


30:58:05.2 


0.0 

-20.5 

0.0 


0.0 

0.0 

-20.5 



continued on following page 



J. Braine et al.: Molecule and Star Formation in M 33 



continued from previous page 



Name 


RA 


Dec 


(5RA 

(") 


(5Dec 

(") 








20.5 


-0.0 








-0.0 


20.5 


M33.7 


01:34:22.37 


30:57:03.1 


0.0 

-20.5 

-0.0 

20.5 

-0.0 


0.0 

-0.0 

-20.5 

-0.0 

20.5 


M33.8 


01:34:29.94 


30:57:10.1 


-20.5 

20.5 

0.0 

0.0 

0.0 


-0.0 
0.0 

20.5 

-20.5 

0.0 


M33.9 


01:34:03.72 


30:55:01.8 


-20.5 

-20.5 

-20.5 

20.5 

20.5 

20.5 

-0.0 

-0.0 

0.0 


-20.5 

-0.0 

20.5 

-20.5 

0.0 

20.5 

20.5 

-20.5 

0.0 


M33.10 


01:34:16.23 


30:52:08.4 


3.7 
0.0 


10.5 
0.0 


M33.il 


01:33:59.46 


30:49:12.1 


-4.6 
0.0 


-28.6 
0.0 


M33.12 


01:34:09.30 


30:49:06.1 


0.0 


0.0 


M33.18 


01:34:13.31 


31:10:32.0 


0.0 

-20.5 

-0.0 

20.5 

0.0 

-0.0 

-10.5 

0.0 
-0.0 
-0.0 
-0.0 


0.0 
0.0 

-20.5 

0.0 

20.5 

-10.5 

-0.0 

-41.0 

-61.5 

-82.0 

-102.5 


M33.20 


01:34:39.79 


31:08:23.8 


0.0 

30.0 

60.0 

90.0 

120.0 

150.0 

180.0 

210.0 

240.0 


0.0 
-0.0 
0.0 
-0.0 
0.0 
-0.0 
0.0 
0.0 
0.0 


M33.22 


01:34:54.09 


31:13:25.5 


0.0 


0.0 


M33.23 


01:34:45.28 


31:12:57.5 


0.0 

22.6 
45.2 
67.8 
90.4 
-10.5 


0.0 

5.6 

11.2 

16.8 

22.4 
0.0 


M33.24 


01:34:52.39 


31:08:36.0 


0.0 


0.0 


M33.25 


01:34:22.69 


31:11:54.4 


0.0 
-0.0 


0.0 
-10.5 


M33.26 


01:34:40.88 


31:13:40.0 


0.0 


0.0 


M33.27 


01:34:51.85 


31:09:39.6 


0.0 


0.0 


M33.34 


01:34:57.40 


31:05:13.9 


0.0 


0.0 



4.2. From isotopic lines and ratios 

Another approach to estimating cloud masses is to measure iso- 
topic lines assumed to be optically thin. In principle, this allows 
one to "count" the molecules of the optically thin species. In 
practice, this is quite difficult but isotopic line ratios provide a 



means to determine whether the GMCs are intrinsically very dif- 
ferent from Galactic clouds. Eight of the clouds have been ob- 
served in '^CO and all of them were detected in the '^CO(l-O) 
line and all but one in '^CO(2-l). 

Assuming optically thin '^CO(l-O) emission and Local 
Thermodynamic Equilibrium at a temperature Tgx >> Tcmb, the 
remarkably narrow range in '^^'''CO line ratios (Table 5) yields 
A^(H2)//co(i-0) values within a facto r 2 of 10^° cm~^/(K km s"') 
(see formulae 14.40 and 14.41 in iRohlfs & WiTsonI |2004 for 
realistic excitation temperatures (assumed 10 - 30 K). This is 
true for both the '^CO(7 =1-0) and 7 = 2-1 transitions. 
The calculations assumed a '^CO abundance with respect to H2 
of (1.2 X 10^) ' - lower '^CO abundances will yield a higher 
A^(H2)//co(i-0) value. Of course, a fraction of the '^CO is prob- 
ably not optically thin and the departure from LTE will be much 
greater in '^CO than '^CO due to the inefficient radiative excita- 
tion. Furthermore, the fraction of the GMC where '''CO is disso- 
ciated will be much greater than for '^CO. The combination of 
these ingredients likely explains the much lower A^(H2)//co(i-0) 
value obtained in this way. Low '^^'''CO line ratios lead to high 
A^(H2)//co(i-0) ratios. 

The '^CO integrated intensities and '^^'^CO line ratios for 
each transition are given in Table 5. The ratios are of or- 
der 10 for the 7 = 2-1 transition and slightly higher for 
the J - 1-0. These values are typi c al of (large) galax- 
ies at large scales (iRickard & Blitz|[T985t ISage & Isbeliiri99ll: 
iPaglione et aI.ll200lT) . Values for individual cl ouds in the MiUcy 
Way are usually lower, in the range 3-5 (iPolk et al.l 119881: 
iDigel et aH 1 1 9941) but may not include all of the low column 
density material with significant '^CO but little '^CO emission. 
Since M33 has a metallicity which is subsolar by a factor 2, the 
'-^/'^CO line ratio could be expected to be higher, as the opti- 
cally thick '^CO will be less affected than the '-^CO by the lower 
abundance. 

The '-^CO spectra shown in Fig. 6 are multiplied by 10 and 
can be seen to lay well over the '^CO spectra. The line ratios 
observed in M33 are thus in good agreement with the picture 
in which molecular clouds in M33 are similar to their Galactic 
counterparts but with a lower metallicity and a A^(H2)//co(i-0) 
higher by about a factor two than in ffie Milky Way disk. 



5. Conclusions 

We present the first detections of both '^CO and '^^CO in the 
extreme outer disk (R > R25) of a subsolar metallicity galaxy. 
Several clouds have been detected in our search for molecular 
gas in the outer disk, enabling us to try to identify the environ- 
ments in which detectable quantities of H2 are found in the dim 
outer disks of spirals. The HI antenna temperature and column 
density and the level of star formation as traced by FIR or Ha 
emission are indicators of the likelihood of detecting CO but 
quite imperfect because two of the strongest outer disk clouds 
would not have been found that way. In and of itself, the unreli- 
ability of this link is interesting. 

As we are close to spatially resolving clouds, the "Virial" 
method is used to estimate molecular cloud masses and 
thus the N{li2)/Ico ratio. The inner disk clouds measured 
are compatible with a Galactic A^(H2)//co factor Our esti- 
mate of A^(H2)//co(i-0) in the outer disk of M 33 is 5 x 
10^" cm"- /(K km s"') with an uncertainty of at most a factor 
2. The '^^'■'CO line ratios are compatible with a A^(H2)//co(i-{)) 
factor somewhat greater than in the Milky Way molecular ring, 
such as the A^(H2)//co ~ 5 x 10^"cm"^/(Kkm s"') estimated 
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Table 2. CO(l-O) detections. 



Table 3. CO(2-l) detectiotis. 



Name 



(5RA 

(") 



(5Dec 

(") 



IcO(l-O) 

(Kkms-') 



rms 
(mK) 



Name 



(5RA 

(") 



(5Dec 



IcO(2-l) 

(Kkms-') 



rms 
(mK) 



M33.2 


0.0 


0.0 


0.20 ± 0.03 


M33.3 


0.0 


0.0 


0.42 ± 0.06 


M33.4 


-10.5 


-10.5 


1.20 ± 0.09 


M33.4 


-10.5 


0.0 


1.50 ±0.08 


M33.4 


-10.5 


10.5 


0.99 ±0.13 


M33.4 


0.0 


-10.5 


1.20 ±0.06 


M33.4 


0.0 


0.0 


1.83 ± 0.05 


M33.4 


0.0 


10.5 


1.25 ±0.08 


M33.4 


10.5 


-10.5 


0.64 ± 0.06 


M33.4 


10.5 


0.0 


0.79 ± 0.06 


M33.4 


10.5 


10.5 


0.60 ± 0.08 


M33.4 


-21.0 


0.0 


0.99 ± 0.07 


M33.4 


0.0 


21.0 


0.94 ± 0.06 


M33.5 


0.0 


0.0 


1.47 ± 0.06 


M33.5 


-20.5 


0.0 


0.24 ± 0.05 


M33.5 


20.5 


0.0 


0.48 ± 0.05 


M33.5 


0.0 


20.5 


0.42 ± 0.07 


M33.7 


0.0 


0.0 


1.50 ±0.04 


M33.7 


-20.5 


0.0 


0.36 ± 0.07 


M33.7 


20.5 


0.0 


0.41 ± 0.07 


M33.7 


0.0 


20.5 


1.04 ±0.06 


M33.8 


-20.5 


0.0 


0.67 ± 0.08 


M33.8 


0.0 


20.5 


0.75 ± 0.07 


M33.8 


0.0 


-20.5 


0.35 ± 0.08 


M33.8 


0.0 


0.0 


1.44 ± 0.04 


M33.9 


-20.5 


-20.5 


0.53 ± 0.08 


M33.9 


20.5 


-20.5 


0.88 ±0.10 


M33.9 


20.5 


0.0 


0.59 ± 0.07 


M33.9 


0.0 


20.5 


1.25 ±0.08 


M33.9 


0.0 


-20.5 


0.74 ±0.15 


M33.9 


0.0 


0.0 


2.51 ±0.04 


M33.10 


3.7 


10.5 


3.46 ±0.14 


M33.10 


0.0 


0.0 


3.55 ± 0.08 


M33.il 


-4.6 


-28.6 


9.99 ±0.17 


M33.il 


0.0 


0.0 


9.73 ± 0.08 


M33.12 


0.0 


0.0 


22.78 ±0.10 


M33.18 


0.0 


0.0 


0.52 ± 0.06 


M33.18 


-20.5 


0.0 


0.26 ± 0.06 


M33.18 


0.0 


-10.5 


0.24 ± 0.03 


M33.18 


-10.5 


0.0 


0.61 ± 0.02 


M33.20 


120.0 


0.0 


0.26 ± 0.03 


M33.20 


150.0 


0.0 


0.33 ± 0.04 



7.1 
17.5 
21.1 
19.1 
21.5 
17.7 
12.5 
16.4 
19.3 
16.9 
17.9 
17.3 
15.9 
14.9 
15.5 
16.2 
16.8 
12.5 
17.9 
16.4 
16.5 
18.3 
18.0 
16.5 

8.1 
25.0 
29.6 
20.6 
22.1 
21.8 
11.7 
28.9 
18.4 
33.5 
16.9 
18.9 
19.5 
19.8 

9.8 

7.6 
11.4 
10.2 



by the Virial method presented here. Adopting this conversion 
factor results in a lower star formation efficiency (SFR / MhJ 
than previously found but the SFE remains higher t han the SFE 
in large local universe spirals (e.g. lKennicut3ll998h . 
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Table 4. Results of gaussian fits to the detected clouds in CO(l- 
0). 

Name velocity line width Mvi, Mi^n, 

(kms-') (kms-') (1000 Mq) (1000 Mq) 



M33.2 

M33.9 

M33.il 

M33.12 

M33.3 

M33.4 

M33.5 

M33.7 

M33.8 

M33.10 

M33.18 



-254.1 : 
-256.2 : 

-251.2: 

-248.2 : 
-258.7 : 
-269.8 : 
-263.0 : 
-263.9 : 

-271.6: 

-267.1 : 
-263.3 : 



0.5 
0.1 
0.1 
0.0 
0.4 
0.1 
0.1 
0.1 
0.2 
0.1 
0.1 



5.0: 
5.3: 
8.5: 

11.4: 
5.5: 
4.4: 
6.7: 
4.1: 

10.7: 
7.6: 
2.5: 



1.2 
0.2 
0.1 
0.1 
1.2 
0.2 
0.3 
0.2 
0.4 
0.2 
0.2 



125 ± 60 
140 ±9 
361 ± 10 
643 ± 1 1 
150 ± 66 

95 ±7 
221 ± 23 

82 ±9 
567 ± 40 
284 ± 17 

31 ±4 



6 

81 

312 

730 

20 

58 

55 

48 

44 

105 

19 



Velocities and line widths, with the errors estimated by the gaussian 
fitting algorithm in CLASS are provided in columns 2 & 3. Making 
the assumption that the clo ud radius is equal to the beamsize, virial 
masses are calculated as in iBraine et alj l ll997ah and the errors using 
only the uncertainty in AV. Masses from the CO(l-O) line intensity are 
given in column 5 using the line area from Table 2a and a "standard" 
galactic conversion A'(H2)//co(i-o) = 2x 10^°cm-^/(Kkm s"'). The 
first 4 clouds have gaussian fits made using the 2.6 km s"' channel 
width whereas the latter 7 fits have been made using data at 0.52 km 
s"' channel width. Cloud 8 appears to have two velocity components, 
producing the broad line. Cloud 2 has a large error associated with 
the CO(l-O) observation; the gaussian fit to the C0(2-l) line yields a 
width of 2.2 ± 0.44A:ra.s"' but the C0(2-l) beam is much smaller and 
may not include all the emission. 
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Table 5. Table 5. '^CO observations 



Name 


1(1-0) 


O"l-0 


1(2-1) 


Hl3J(l-0) 


I(if)(2-1) 




Kkm.r' 


mK 


Kkmi--' 






M33.4 


0.15 ±0.02 


5.6 


31.2 mK 


12.2 





M33.7 


0.14 + 0.02 


4.2 


0.16 ±0.01 


10.7 


8.5 


M33.8 


0.13 ±0.02 


3.3 


0.23 ± 0.02 


11.2 


10.0 


M33.9 


0.16 ±0.01 


3.9 


0.23 ± 0.04 


15.7 


10.1 


M33.10 


0.29 ± 0.02 


5.2 


0.22 ± 0.02 


12.2 


13.8 


M33.il 


0.77 ± 0.02 


5.1 
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12.6 


9.8 


M33.12 
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4.6 
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8.9 


11.4 
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Fig. 2. CO(l-O) spectra (yellow) and HI spectra (red dotted) of sources M33_20, M33_24, and M33_27 on 70//m image with HI 
contours at 0.8, 1.2 (dashed), and 1.6 xlO^' cm"^ (dotted). The wedge shows the VO/zm image scale in MJy sr"'. Note that two 
positions along the strip (M33_20) have been detected but that they correspond to neither HI nor IQjim maxima. 
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Fig. 3. CO(l-O) spectra (yellow) and HI spectra (red dotted) of sources M33_22, M33_23, and M33_26 on 70//ot image with HI 
contours at 0.8, 1.2 (dashed), and 1.6 xlO^' cm"^ (dotted). The wedge shows the IQ^m image scale in MJy sr"'. None of the 
positions have been detected in CO. 



